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The superconducting state of an infinitely long supercon- 
ducting cylinder surrounded by a medium which enhances its 
superconductivity near the boundary is studied within the 
nonlinear Ginzburg-Landau theory. This enhancement can 
be due to the proximity of another superconductor or due to 
surface treatment. Quantities like the free energy, the mag- 
netization and the Cooper-pair density are calculated. Phase 
diagrams are obtained to investigate how the critical field 
and the critical temperature depend on this surface enhance- 
ment for different values of the Ginzburg-Landau parameter 
k. Increasing the superconductivity near the surface leads 
to higher critical fields and critical temperatures. For small 
cylinder diameters only giant vortex states nucleate, while for 
larger cylinders multivortices can nucleate. The stability of 
these multivortex states also depends on the surface enhance- 
ment. For type-I superconductors we found the remarkable 
result that for a range of values of the surface extrapolation 
length the superconductor can transit from the Meissner state 
into superconducting states with vorticity L > 1. Such a 
behaviour is not found for the case of large ft, i.e. type-II 
superconduct ivity. 
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I. INTRODUCTION 

The study of the superconducting properties of meso- 
scopic samples received recently a lot of attention due to 
the progress in nanofabrication technologies. A sample 
is mesoscopic if its size is comparable to the penetration 
depth A and/or the coherence length £. 

The properties of superconductors and therefore their 
applications are determined by the critical parameters, 
i.e. the critical fields (H c , H c \, H C 2, H C 3), the critical 
current j c and the critical temperature T c . For applica- 
tion purposes one wants to improve these parameters. In- 
stead of searching for new bulk materials with new phys- 
ical properties, one can also try to modify the properties 
of an existing material by nanostructuring the supercon- 
ductor, i.e. 'quantum design'. Therefore, theoretical and 
experimental studies directed to improve the critical pa- 
rameters of mesoscopic superconductors are important. 

From the point of view of improving these parameters, 



mesoscopic thin disks immersed in an insulating medium 
and placed in a perpendicular magnetic field have at- 
tracted most attention. In disks with a small radius only 
giant vortex states exist which are axially symmetric su- 
perconducting states Jl)-^. The different giant vortex 
states are characterized by their angular momentum, i.e. 
their vorticity L. For disks with a sufficiently large ra- 
dius, multivortex states can exist which are the equiva- 
lent of the Abrikosov vortex lattice in bulk type-II su- 
perconductors (k = A/£ > l/v2) and the vorticity cor- 
responds to the number of separate vortices [p|-p8[. For 
such disks, it is found that the critical magnetic field de- 
pends strongly on the size and on the Ginzburg-Landau 
parameter and less on the thickness of the disk. 

Fink and Prcsson |H]j2(| studied the Meissner state 
and the giant vortex state of a superconducting infi- 
nite cylinder immersed in an insulator or vacuum and 
placed in a magnetic field parallel to the cylinder axis 
using the Ginzburg-Landau theory. They calculated how 
the critical field depends on the cylinder radius and 
on the Ginzburg-Landau parameter k. Zharkov et al 
Ipif studied similar superconducting cylinders using the 
Ginzburg-Landau theory and calculated the critical fields 
as a function of the cylinder radius, the temperature and 
the Ginzburg-Landau parameter. 

Varying the sample size and the Ginzburg-Landau pa- 
rameter, one can decrease or increase the critical field, 
and thus the critical current, but it does not influence 
the critical temperature. How can one increase the crit- 
ical temperature? Fink and Joiner considered a 
semi-infinite superconducting half-space where the sur- 
face was treated by cold working in such a way that 
the superconductivity near the surface was enhanced, 
i.e. the slope of the superconducting order parameter 
increased near the sample surface. They found that such 
a surface treatment leads to higher critical temperatures, 
larger critical fields, and larger critical currents. Another 
way to enhance superconductivity near the surface is by 
bringing the superconductor in contact with a well cho- 
sen superconducting or semiconducting layer |23| . The 
added superconducting layer must have a higher transi- 
tion temperature than the superconducting sample and 
in this case the surface enhancement of superconductiv- 
ity is caused by the proximity effect. It is necessary that 
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the semiconducting layer has a bandgap which overlaps 
the superconducting gap of the superconductor. 

Montevecchi and Indekeu |24| performed a theoretical 
study of the effect of confinement on the superconduct- 
ing/normal transition for systems with surface enhance- 
ment. Using the linear Ginzburg-Landau theory they 
studied the critical temperature at zero field for a thin 
film, an infinite cylinder and a sphere. For all geome- 
tries they found that at H = the critical temperature 
increases with the enhancement of surface superconduc- 
tivity. For thin films in a parallel field, they also calcu- 
lated a H — T phase diagram. They found that surface 
enhancement leads to an increase of the critical field H c s 
and of the critical temperature. Notice that their results 
do not depend on k, since they use the linear theory. 

Yampolskii and Peeters |2f| investigated theoretically 
the vortex structure of thin mesoscopic disks in a perpen- 
dicular magnetic field surrounded by a medium which en- 
hances surface superconductivity. If the size of these thin 
disks is sufficiently large, they found giant vortex states 
as well as multivortex states appearing as metastable 
states and as ground states. The enhancement of the 
superconductivity near the surface leads to a stabiliza- 
tion of the multivortex states as ground states. They 
also calculated a H — T phase diagram which showed 
that the critical temperature and the critical field were 
significantly increased by enhancing the surface super- 
conductivity. 

In the present paper we investigate the effect of the 
enhancement of surface superconductivity on the criti- 
cal field and the critical temperature for superconducting 
cylinders with radii equal to a few coherence lengths £. 
We also study the influence of the Ginzburg-Landau pa- 
rameter k. Our theoretical analysis is based on a full self- 
consistent numerical solution of the coupled non-linear 
Ginzburg-Landau equations. No a priori arrangement of 
the vortex configuration or of the type of vortex config- 
uration is assumed. 

The paper is organized as follows. In Sec. II we present 
our theoretical model. We explain how the enhancement 
of superconductivity near the surface is taken into ac- 
count. In Sec. Ill we consider small cylinders, i.e. in- 
finitely long superconducting cylinders with small radius. 
In such small cylinders only axially symmetric states or 
giant vortex states nucleate. We discuss how the criti- 
cal field and the critical temperature are influenced by 
the surface enhancement of superconductivity. In Sec. 
IV larger cylinders are studied where multivortices ap- 
pear. We investigate the dependence of the nucleation of 
these states on the surface enhancement. Our results are 
summarized in Sec. V. 



conductivity at the edge of the cylinder. Along the axis 
of the cylinder a uniform magnetic field Hq is applied. 
To deal with this problem we use the Ginzburg-Landau 
theory and we solve numerically and self-consistently the 
system of two coupled equations 
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where the density of superconducting current j is given 
by 
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The boundary condition for the vector potential is such 
that far away from the cylinder the field equals the ex- 
ternal field Hq, i.e. 



A 



r — >oo 2 



(2a) 



where p indicates the radial position from the center of 
the cylinder and "e^ is the azimuthal direction. The 
general boundary condition for the order parameter can 
be written as 
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where ft is the unit vector normal to the cylinder and b is 
the surface extrapolation length. The value of b is deter- 
mined by the medium in which the cylinder is immersed. 
If the surrounding medium is vacuum or an insulator one 
has b — > oo, for metals b > and for ferromagnets 6^0. 
When, for example, the surface of the cylinder is in con- 
tact with a superconducting layer of a higher T c one has 
b < 0. In the present paper we will consider the latter 
case of b < 0. 

Using dimensionless variables and the London gauge 
divA = the GL equations can be written as 
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and the boundary condition for the order parameter as 



II. THEORETICAL FORMALISM 

We consider infinite superconducting cylinders with ra- 
dius R surrounded by a medium which enhances super- 
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Here the distance is measured in units of the coher- 
ence length £, the order parameter in "to = \/—ot//3, 
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the vector potential in ch/2e£, and the magnetic field in 
H c2 = ch/2e( 2 = kV2H c , where H c = v /4vra 2 //3 is the 
critical field. 



To solve the system of Eqs. (^a-3b), we apply a finite- 
difference representation of the order parameter and the 
vector potential on a uniform Cartesian space grid (x,y), 
with typically 128 grid points over a distance 2R Q (i.e. 
the diameter of the cylinder), and we used the link vari- 
able approach |27|], and an iteration procedure based on 
the Gauss-Seidel technique to find \P. The vector po- 
tential is obtained with the fast Fourier transform tech- 
nique where we set A \ x \ = r s y \ y \ = n s = Hq (x, —y) /2 at the 
boundary of a larger space grid (typically Rs — 4i?„). 

Two different vortex configurations are possible in ax- 
ial symmetric samples; the giant vortex state and the 
multivortex state. The giant vortex state has cylin- 
drical symmetry, which leads to the order parameter 
$ = ip (p) exp{iLcf>), where p and cj) are the cylindrical 
coordinates and L is the angular momentum or vorticity. 
Along a closed path, which lies near the edge of the cylin- 
der, the order parameter phase difference is given by L 
times 2-7T. A multivortex state is described by a mixture 
of different angular harmonics which is possible due to 
the nonlinearity of the GL equations. In this case the 
vorticity is nothing else than the number of vortices in 
the superconductor. 

To find the different vortex configurations, which in- 
clude the metastable states, we search for the steady- 
state solutions of Eqs. j3a[-^h|) starting from different 
randomly generated initial conditions. Then we in- 
crease / decrease slowly the magnetic field and recalculate 
each time the exact vortex structure. We do this for each 
vortex configuration in a magnetic field range where the 
vorticity stays the same. 

In order to determine which vortex state corresponds 
to the ground state we calculated the difference in free 
energy density between the superconducting state and 
the normal state 
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where I? = rot A, V is the volume of a cylinder with 
unit height, S the surface of the infinite cylinder and 
F = The contribution of the surface is taken 

into account by the last term of Eq. (0) . By comparing 
the dimensionless free energy of the different giant and 
multivortex configurations, we obtain the ground states 
and the metastable states. The dimensionless magneti- 
zation is defined as 
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where (H) is the magnetic field averaged over the sample. 
The magnetization is a direct measure of the magnetic 
field expelled from the sample. 



The temperature is included in £, A, H c2 , through their 
temperature dependencies 
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where T c q is the critical temperature at zero magnetic 
field for the normal boundary condition, i.e. — £(0)/6 = 
0. We will only insert the temperature explicitly if we 
consider the H—T phase diagrams, while the other calcu- 
lations are for (arbitrary) fixed temperature. This means 
that we vary b at constant £ and constant R. Notice fur- 
ther that the Ginzburg-Landau parameter n — A/£ is 
independent of the temperature. 



III. SMALL CYLINDERS 

First we discuss infinitely long cylinders with a small 
radius, where we may limit ourselves to the giant vortex 
states. In this case the confinement effects are dominant 
and they impose the cylindrical symmetry of the bound- 
ary of the cylinder on the vortex state. Consequently, 
the dimensions of the Ginzburg-Landau equations are re- 
duced, which improves the accuracy and the computation 
time. 



A. Type-I cylinders 

We consider infinite long cylinders with radius R = 
2.0£ for a Ginzburg-Landau parameter k — 0.28. 

Figs. [j] (a,b) show the ground state free energy, i.e. the 
thermodynamic equilibrium free energy, and the magne- 
tization as a function of the applied magnetic field for 
such a superconducting cylinder for the usual bound- 
ary condition, i.e. — £/b = 0.0 (solid curves), and for 
— £/6 = 0.2 (dashed curves) and 0.4 (dash-dotted curves), 
which correspond to surface enhancement. With in- 
creasing — £/& the free energy at zero magnetic field be- 
comes more negative, which means an enhancement of 
superconductivity. Also the ground state superconduct- 
ing/normal transition moves to higher fields. Notice that 
the magnetic field expulsion (see Fig. 0(b)) is enhanced 
with increasing — £/&. 

Figs. ||(a-c) show the radial dependence of the Cooper- 
pair density, the magnetic field and the current density, 
respectively, for the previously considered cylinder at the 
applied magnetic field H /H c2 — 2.02 and for vorticity 
L = 0. The solid curves give the result for — £/6 = 0, 
the dashed curves for — £/& = 0.2 and the dash-dotted 
curves for — £/6 = 0.4. From Fig. ||(a) it is very clear 
that the Cooper-pair density near the surface, and hence 
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the surface superconductivity, increases with increasing 
— £/6. For a small cylinder radius, also the Cooper-pair 
density in the center of the cylinder is influenced by the 
boundary condition. From Fig. ^|(b) one can see that en- 
hancing the surface superconductivity results in a more 
pronounced Meissner effect, i.e. the magnetic field ex- 
pulsion becomes more complete with increasing — £/&. 
Notice that the magnetic field for this situation is al- 
ways zero in the center of the cylinder. To expel the 
magnetic field more, the superconductor has to induce 
more superconducting current near the surface. Fig. ^|(c) 
shows that the current density is zero in the center of 
the superconductor and becomes more negative near the 
surface. The superconducting Meissner currents are also 
more concentrated near the surface with increasing — £/&. 
For the Meissner state the current has the same sign in 
the whole sample and is thus flowing in the whole sample 
in the same direction. 

In Fig. |H the dependence of the transition fields for 
the different giant vortex states on the surface enhance- 
ment b is given for a superconducting cylinder with ra- 
dius R — 2.0£ and k = 0.28. The solid curves give the 
ground state transitions between the different L-states 
and the thick solid curve indicates the superconduct- 
ing/normal transition. For the normal boundary con- 
dition, — £/6 = 0, it is evident that the ground state 
consists of the Meissner state, i.e. L = 0, since K is suf- 
ficiently small. The superconducting/normal transition 
occurs at Hq/H C 2 ~ 2.9. What happens with increas- 
ing — £/&? (i) For — £/b < 0.65 the ground state is still 
given by the Meissner state, and with increasing — £/6 
the superconducting/normal-transition moves to higher 
fields, (ii) For — £/b > 0.65 different L states become 
the ground state. This is in agreement with the results 
of Montevecchi |28j who also found transitions between 
different L states for 1/6 < for a fixed cylinder radius 
Rj \b\ = 1.2 and re = 0.3 (see Fig. 1.14 of Ref. @). Re- 
markably we find that the ground state does not evolve 
from the Meissner state to a state with vorticity L = 1. 
Rather it evolves to a state with larger vorticity. For 
0.65 < — £/6 < 0.8 the ground state changes from the 
Meissner state to a state with L = 4, for 0.8 < — £/b < 1.8 
to a state with L — 3, and for — £/b > 1.8 to a state with 
L = 2. (iii) Further increasing — £/6 the ground state will 
change first from L = to L = 1. The latter transition 
corresponds to type-II behavior which only occurs in the 
absence of surface enhancement for k larger than some 
critical Ginzburg-Landau parameter k,i . 

The three different regimes which we found for the case 
of increasing — £/b for fixed small n are very similar to the 
three regimes for the case of increasing n for the normal 
boundary condition, i.e. — £/6 = 0: (i) For the normal 
boundary condition it is known that only the Meiss- 
ner state is the ground state for k smaller than a critical 
parameter k\. (ii) For n\ < k < k% surface supercon- 
ductivity can occur and states with vorticity larger than 
one exist. For example, Fink and Presson Jl9| found 
ground state transitions from the Meissner state with 



L = into surface superconducting states with L = 4 
for a superconducting cylinder with radius R = 3.0£ and 
k = 0.5. (iii) For k > k 2 (type II superconductivity) the 
ground state transits from the Meissner state to states 
with L = 1, 2, n, successively, where the value of n 
depends on the radius and the Ginzburg-Landau param- 
eter (see for example Ref. |l^]). In bulk superconductors 
Ki w 0.42 and n-i = l/y2, but for cylinders these param- 
eters are radius dependent Q . 

What is the reason for this remarkable behaviour of the 
ground state for increasing — £/&? (i) For — £/6 < 0.65 the 
order parameter, and thus the Cooper-pair density, near 
the cylinder boundary increases with increasing — £/b 
which also leads to an enhancement of the order pa- 
rameter in the center, albeit less pronounced (see also 
Fig. ||(a)). Consequently, the free energy at zero mag- 
netic field becomes more negative and the superconduct- 
ing/normal transition moves to higher fields (see also 
Fig. (a)), (ii) With increasing — £/b the Cooper-pair 
density in the center for the L = state is not increasing 
as strongly as the Cooper-pair density near the bound- 
ary and the cylinder remains superconducting at higher 
fields. For — £/b > 0.65 it becomes energetically less 
favourable to expell rather high fluxes from the sam- 
ple. The effect of the increase of vorticity on the free 
energy is shown in Fig. |4| for — = 0.7. From the in- 
set it is clear that the free energy of the L = 4 state 
becomes lower than the one of the Meissner state at 
Hq/H C 2 Rj 5.12. This means that for this field the pene- 
tration of 4 fluxes becomes more favourable than the ex- 
pulsion of the field. At H /H c2 « 5.20, 5.30 and 5.32 the 
ground state changes respectively into the L = 5 state, 
the L = 6 state and the normal state. Further increasing 
— £/6 leads to higher Cooper-pair density near the bound- 
ary and, consequently, to higher induced supercurrents. 
As a result, a smaller amount of flux will penetrate into 
the sample after the first ground state transition. For 
example, with increasing field at — £/6 > 0.8 the ground 
state changes from the Meissner state into the L = 3 
state and at higher fields into states with L = 4, 5, 6 
and so on. (iii) For high values of — £/b, the explanation 
is analogous as above but now the ground state changes 
from L = to L = 1, 2, 3 and so forth. 

In Fig. H we plot the H — T phase diagram for a su- 
perconducting cylinder with R = 2.0£ and K = 0.28 
for -£(0)/& = 0.0, 0.1 and 0.2. For these values of 
— £(0)/& the superconducting ground state is always given 
by the Meissner state (see also Fig. |). The supercon- 
ducting/normal transition is shown by solid curves for 
-f(0)/6 = 0.0, by dashed curves for -£(0)/6 = 0.1, and 
by dash-dotted curves for — £(0)/6 = 0.2. For k = the 
superconducting/normal transition is a straight line in 
the H — T phase diagram. For k — 0.28, the supercon- 
ducting/normal transition is still given by a straight line 
for T <C T c , regardless of the value of — £(0)/6, but it 
has now a curvature near T c . With increasing — £(0)/6 
the superconducting/normal transition, which is of first 
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order, moves to higher temperatures for fixed field or to 
higher fields for fixed temperature. 

B. Type-II cylinders 

We consider infinite cylinders with radius R = 2.0£ for 
a Ginzburg-Landau parameter k = 1.0. 

Figs. ^(a,b) show the ground state free energy and mag- 
netization as a function of the applied magnetic field 
for such a superconducting cylinder for — £/b = 0.0 
(solid curves), 0.2 (dashed curves) and 0.4 (dash-dotted 
curves). Now, the superconducting ground state is not 
always the Meissner state. For example, for the usual 
boundary condition, i.e. — £/& = 0.0, the vorticity of 
the ground state is L = for H /H C 2 < 1.195. At 
the first transition field Hq/H C 2 = 1.195 the vorticity 
of the ground state changes from L = to L = 1 and 
then it remains L = 1 until H /H C 2 = 1.7575 where it 
changes to L = 2. At Hq/H C 2 = 2.12 the free energy 
of the ground state becomes zero and the superconduct- 
ing/normal transition takes place. The transitions be- 
tween the different L-states are indicated by corners in 
the free energy and jumps in the magnetization. As for 
type-I cylinders, the free energy at Hq — becomes more 
negative with increasing — £/&, i.e. superconductivity is 
enhanced, and the transition to the normal state moves 
to higher fields. Furthermore, the vorticity of the ground 
state can become larger and the peaks in the magnetiza- 
tion are higher, indicating a more efficient expulsion of 
the magnetic field from the superconductor. Notice that 
the magnetic fields for the transitions between different 
L-states are almost independent of the value of —£/b, 
but the surface critical field H c3 is a sensitive function of 

Next, we investigate the radial dependence of the 
Cooper-pair density, the magnetic field and the current 
density for such a superconducting cylinder. For vor- 
ticity L = the results are analogous to the results 
for type-I cylinders which were described in Fig. 0. In 
Figs. |7](a-c) we show the radial dependence of, respec- 
tively, the Cooper-pair density, the magnetic field and the 
current density for the ground state at H^/H C 2 = 2.02, 
i.e. the giant vortex state with vorticity L = 2. The 
solid curves give the result for — £/b = 0, the dashed 
curves for — £/b = 0.2 and the dash-dotted curves for 
— £/6 = 0.4. In Fig. 0(a) it is shown that the Cooper- 
pair density is zero for all — £/b in the center of the cylin- 
der, i.e. at the position of the giant vortex with vorticity 
L = 2. Near the surface the Cooper-pair density is en- 
hanced, which is understandable because the slope of 
at p = R is given by —£/b. For the giant vortex state, the 
magnetic field is not zero in the center, as was the case 
for L = 0, but it can be even higher than the external 
field Hq = 2.02iJ C 2. With increasing — £/&, the supercon- 
ductivity near the boundary is enhanced and, thus, more 
magnetic field can be expelled from the cylinder and the 



giant vortex will be more compressed in the center (see 
Fig. 0(b)). As a consequence, the minimum of magnetic 
field inside the cylinder decreases and the magnetic field 
at the position of the vortex increases. To expel the mag- 
netic field more and to compress the giant vortex better 
in the center, the superconductor has to induce a larger 
superconducting current. Therefore the superconducting 
current is more positive close to the giant vortex and 
more negative near the boundary (see Fig. 0(c)). 

The fact that the transition fields between different in- 
states are almost independent of the value of —£/£>, can 
also be seen from the (—£/&) — H phase diagram for a su- 
perconductor cylinder with radius R = 2.0£ and K = 1.0, 
which is shown in Fig. The thin solid curves give the 
ground state transitions between the different L-states 
and the thick solid curve indicates the superconduct- 
ing/normal transition. For k = 1.0 we do not find ground 
state transitions between a state with vorticity L — and 
vorticity L > 1, as was the case for n = 0.28. Moreover, 
the magnetic field range over which the ground state has 
a particular vorticity L is almost the same for all L > 0, 
namely AH L = H L ^ L+1 - E L -x^L « 0.57iJ c2 . Fig. 
indicates that the width of the superconducting sheath 
in p/£ is independent of b. The value of b influences 
only the amplitude of the order parameter near the sur- 
face (Fig. @(a)), the diamagnetization (Fig. 0(b)) and the 
strength of the currents (Fig. 0(c)). Therefore, the space 
available for the confined flux in the core region of the 
cylinder does (almost) not depend on the enhancement. 
As a consequence, the number of fluxoids, or vorticity Z, 
depends only on the strength of the external field Ho and 
almost not on b for fixed temperature. 

In Fig. we plot the H — T phase diagram for a su- 
perconducting cylinder with R = 2.0£ and k = 1.0 for 
— £(0)/6 = 0.0, 0.1 and 0.2. The ground state transitions 
are shown by solid curves for — £(0)/6 = 0.0, by dashed 
curves for — £(0)/& = 0.1, and by dash-dotted curves for 
— £(0)/6 = 0.2. The thick curves indicate the supercon- 
ducting/normal transition and the thinner curves (they 
are almost straight lines) show the ground state transi- 
tions between the giant vortex states L L + 1. From 
Fig. it is clear that the transitions between the different 
L states are (almost) independent of the value of ~£,/b, 
but the superconducting/normal transition is a sensi- 
tive function of — £/6 (see also Fig. 0). With increasing 
— £(0)/6 the superconducting/normal transition moves to 
higher temperatures for fixed field or to higher fields for 
fixed temperature as was also the case for n = 0.28. The 
zero field critical temperature is the same for k = 1.0 as 
for k — 0.28 and depends only on — £(0)/6. The criti- 
cal fields at zero temperature are different in both cases. 
For -£(0)/6 = 0.2, for example, H c3 (0) = 3.5£T c2 (0) 
for k = 0.28 and H c3 (0) = 2.7i/ c2 (0) for k = 1.0. In 
fact, the 'zero' temperature corresponds to some refer- 
ence temperature rather than the real temperature be- 
cause the Ginzburg-Landau theory is not valid at very 
low temperatures. Another significant difference is that 
for k — 1.0 the superconducting/normal transition is 
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no longer a straight line for temperatures T <C T c , but 
consists of corners which indicate the transition between 
the different L states. These corners are bicritical points 
where three phases coexist for a given value of — £(0)/6. 
For example, the normal state coexists with the L = 
state and the L = 1 state at Ho/H c2 (0) = 0.86 and 
T/T c0 = 0.84 for -f(0)/& = 0.2. Note that all tran- 
sitions between different L states are of first order and 
correspond to free energy crossings, while the supercon- 
ducting/normal transition is a second order transition. 
These bicritical points occur at magnetic fields which are 
only weakly sensitive to the surface enhancement, while 
the transition temperatures are quite sensitive to b. No- 
tice further that the superconducting/normal phase tran- 
sition boundary is exactly the same for a cylinder with 
k = 1.0 as for a very thin disk (thickness d -C £) with the 
same radius, if the same boundary condition (Eq. (2b)) 
is taken on the sides, but the standard boundary condi- 
tion (— £/& = 0) on the top and the bottom of the disk. 
This leads to a z-independent order parameter and, as a 
consequence, to similar physics as for an infinitely long 
cylinder. This can be seen by comparing Fig. ^| with 
Fig. 14 of Ref. In the case of the disk, the field is 

perpendicular to the disk. 



IV. LARGE CYLINDERS: MULTIVORTICES 

Up to now, we considered only cylinders with small 
radii where the confinement effects were dominant and 
only giant vortex states were stable. Now, we consider su- 
perconducting cylinders with a large radius in which mul- 
tivortex states can nucleate for certain magnetic fields. 
This means that we can no longer limit our calculations 
to axially symmetric solutions. Nevertheless, from the 
study of disks [|| and rings jL4| we know that the transi- 
tion fields between states with different vorticity do (al- 
most) not depend on the fact whether one considers axial 
symmetry or not. Therefore, we can still calculate the 
(—£/&) — H phase diagrams using the cylindrical sym- 
metry approach, but if we want to know the real vortex 
configurations, we have to use the general theory in order 
to account for multivortices. 

We calculate the (—£/&) — H phase diagrams for a su- 
perconducting cylinder with radius R = 4.0£ for k = 0.28 
and K — 1.0, which are shown in Figs. [l0|(a,b), respec- 
tively. The thin solid curves give the ground state tran- 
sitions between the different L-states and the thick solid 
curve indicates the superconducting/normal transition. 
The conclusions are similar to those for small cylinders, 
but with many more transitions between different in- 
states. However, for k = 0.28 we find now transitions 
from the Meissner state immediately to a state with vor- 
ticity L — 13 as a function of the magnetic field in the 
range 0.43 < — £/b < 0.47. With further increasing — £/b 
we find first transitions from L = to L = 12, 11, 10 and 
9, respectively. Notice that triple points occur where the 



Meissner state coexists with two different giant vortex 
states. For example, at Hq/H C 2 ~ 3.4 and — £/6 w 0.47 
the Meissner state coexists with the L = 12 state and the 
L = 13 state. This means that these three states have 
the same free energy in that case. 

For k = 1.0 the magnetic field range over which the 
ground state has a particular vorticity L is almost the 
same for all L > 0, namely AHl ss 0.125ii C 2, which 
means that the ground state changes more quickly from 
vorticity L than for cylinders with radius R — 2.0£, where 
AHl ~ 0.57 H c2 , which is due to the larger cross section 
of the cylinder. 

Next we will compare the free energy and the mag- 
netization for the usual boundary condition — £/b = 
with the case of surface enhancement, — £/b = 0.2. 
Figs. ll(a,b) show the free energy and the magnetiza- 
tion for a cylinder with radius R — 4.0£ and k = 1.0 
for the usual boundary condition. The different giant 
vortex states are indicated by the solid curves, the mul- 
tivortex states by the dashed curves and the multivortex 
to giant vortex transition fields Hug by the open cir- 
cles. Notice that we find superconducting states up to 
vorticity L = 11 and that the superconducting/normal 
transition field is H c3 /H c2 ~ 1.85 (see also Fig. |lO|(b)). 
Multivortices can nucleate for vorticities L = 2, 3, and 
4. The multivortex state with L = 2 can nucleate for 
0.295 < H /H c2 < 0.92, with L = 3 for 0.42 < H /H c2 < 
1.045, and with L = 4 for 0.5325 < H /H c2 < 0.92. In 
Fig. [ll](b) the equilibrium ground state transitions are 
indicated by the vertical curves. 

Figs. |l^(a,b) show the free energy and the magnetiza- 
tion for a cylinder with radius R = 4.0£ and k — 1.0 for 
— £/6 = 0.2. The different giant vortex states are indi- 
cated by the solid curves, the multivortex states by the 
dashed curves and the multivortex to giant vortex transi- 
tion fields Hmg by the open circles. The latter transition 
is of second order. For this boundary condition the su- 
perconducting/normal transition field is H c ^,/H c2 « 2.35 
which is appreciably higher than for — £/b = 0. Also 
vortex states with higher vorticity are stable, up to 15 
instead of up to 11 for — £/& = 0. Multivortex states can 
nucleate with L = 2 for 0.1325 < H Q /H c2 < 0.945, with 
L = 3 for 0.2575 < H /H c2 < 1.02, and with L = 4 for 
0.3825 < H Q /H c2 < 0.895. 

How does the value of — £/b influence the nucleation 
of the multivortex states? Fig. |l3| shows the free energy 
of the L — 3 state for a superconducting cylinder with 
radius R = 4.0£ and k = 1.0 for -f/6 = 0.0, 0.2, 0.4, 0.6, 
0.8 and 1.0 as a function of the applied magnetic field. 
Giant vortex states are given by solid curves and multi- 
vortex states by dashed curves. The transitions from a 
multivortex state to a giant vortex state are indicated by 
the open circles. With increasing — £/b the L = 3 state 
can nucleate up to larger fields, i.e. H nuc /H c2 — 1.21, 
1.32, 1.47, 1.66, 1.88 and 2.16 for -£/b = 0.0, 0.2, 0.4, 
0.6, 0.8 and 1.0, respectively. On the other hand the 
multivortex to giant vortex transition field decreases, i.e. 
H M g/H c2 = 1.045, 1.02, 1.0075, 0.995, 0.97 and 0.82 for 
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-i/b = 0.0, 0.2, 0.4, 0.6, 0.8 and 1.0, respectively. This 
means that with increasing — £/6 the magnetic field re- 
gion over which the giant vortex state exists increases and 
the one of the multivortex state decreases. The reason 
is that increasing — £/& corresponds to an enhancement 
of superconductivity near the boundary. Therefore, it is 
more difficult for the magnetic field to penetrate near the 
cylinder edge than in the center. Since multivortices are 
situated more closely to the cylinder edge, it is obvious 
that multivortices become unstable at lower fields for in- 
creasing —£/b. Or alternatively, the confinement effects 
become more important for increasing surface supercon- 
ductivity. 

To illustrate this better, we show the Cooper-pair den- 
sity of the L = 3 state for such a cylinder in Figs. |l4|(a-f) 
at H /H c2 = 0.445 for -f/6 = 0.0, 0.2, 0.4, 0.6, 0.8 and 
1.0, respectively. High Cooper-pair density is given by 
dark regions, low Cooper-pair density by light regions. 
The three vortices correspond to the three white spots. 
For the usual boundary condition — £/b = 0.0 the vortices 
are clearly separated. With increasing — £/b the vortices 
are pushed to the center and at — £/b = 0.4 they already 
start to overlap. For — £/b = 0.8 and 10 the vortices are 
very close to each other. Therefore, we plotted the log- 
arithm of the Cooper-pair density to show the positions 
of the vortices and to prove that it it still a multivortex 
state. Thus, by the enhancement of the Cooper-pair den- 
sity near the boundary the vortices move to the center 
and for increasing — £/6 they will recombine in the center 
creating a giant vortex state. 

V. CONCLUSIONS 

In the present paper we investigated the effect of the 
enhancement of surface superconductivity on the criti- 
cal field and the critical temperature for superconducting 
cylinders with radii comparable to the coherence length £. 
We also studied the influence of a the Ginzburg-Landau 
parameter k. A distinction was made between cylinders 
with small radii where the confinement effects dominate 
and only giant vortex states exist, and cylinders with a 
larger radius where multivortices can nucleate for certain 
magnetic fields and vorticitics. 

Generally, increasing — £/b leads to a more negative 
free energy at H = and to a higher superconduct- 
ing/normal transition field. We also studied the mag- 
netic field distribution, the Cooper-pair density and the 
current density for different values of the surface enhance- 
ment. For higher values of — £/b more magnetic field can 
be expelled from the cylinder and the giant vortex in 
the center is compressed more. Therefore, higher cur- 
rents are induced near the boundary and near the giant 
vortex. The Cooper-pair density close to the boundary 
increases, and for small cylinders this also influences the 
Cooper-pair density in the center. 

From the H — T phase diagrams we found that the 



critical temperature depends on —£/b, while it is inde- 
pendent of k (cf. Ref. pij] ). With increasing — £/b the 
critical temperature increases for fixed magnetic field and 
the critical magnetic field increases for fixed temperature. 
We also obtained (—£/&) — H phase diagrams. In type- 

I cylinders the surface enhancement has drastic conse- 
quences. Even at low k the surface enhancement leads to 
transitions between different L-state and thus to type- 

II behavior. Moreover, as a function of the magnetic 
field the superconducting ground state transits from the 
Meissner state to a vortex state with L > 1 over a range 
of — £/b values. Therefore, we can conclude that increas- 
ing — £/b for type-I superconductors seems to have a simi- 
lar effect as increasing k for fixed — £/b for certain proper- 
ties of the superconducting cylinder. In type-II cylinders 
we found that the magnetic field range over which the 
ground state has a particular vorticity L is almost the 
same for all L > 1. This magnetic field range decreases 
with increasing cylinder radius. 

If the cylinder radius is sufficiently large multivortex 
states can nucleate and we studied the influence of the 
surface enhancement on the nucleation of these states. 
We found that at fixed field and with increasing — £/b the 
multivortices move to the center creating a giant vortex 
state. Thus surface enhancement destabilizes the multi- 
vortex state. 
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FIG. 2. (a) The Cooper-pair density, (b) the magnetic field, 
and (c) the current density as a function of the radial position 
for a cylinder with radius R — 2.0£ and for k = 0.28. The 
external field is Hq/H C 2 = 2.02 and the vorticity is L — 0. 

FIG. 3. The (—£/&) — H phase diagram for a superconduct- 
ing cylinder with radius R = 2.0£ and k = 0.28. The solid 
curves indicate the transitions between the different Instates 
and the thick curve gives the superconducting/normal tran- 
sition. 

FIG. 4. The free energy for a superconducting cylinder 
with radius R/£ = 2.0 and k = 0.28 for -£/& = 0.7 as a 
function of the external magnetic field. 

FIG. 5. The H — T phase diagram for a superconducting 
cylinder with radius R = 2.0£ and n = 0.28 for -£/& = 0.0 
(solid curves), 0.2 (dashed curves) and 0.4 (dash-dotted 
curves) . 

FIG. 6. (a) The ground state free energy and (b) the mag- 
netization as a function of the applied magnetic field for a 
superconducting cylinder with radius R = 2.0£ for k = 1.0 
and — £/& = (solid curves), 0.2 (dashed curves) and 0.4 
(dash-dotted curves). 

FIG. 7. (a) The Cooper-pair density, (b) the magnetic field, 
and (c) the current density as a function of the radial position 
for a cylinder with radius R — 2.0£ and for k = 1.0. The 
external field is Hq/H C 2 = 2.02 and the vorticity is L — 2. 

FIG. 8. The (— £/&) — H phase diagram for a supercon- 
ducting cylinder with radius R = 2.0£ and k = 1.0. The solid 
curves indicate the transitions between the different Instates 
and the thick curve gives the superconducting/normal tran- 
sition. 

FIG. 9. The H — T phase diagram for a superconduct- 
ing cylinder with radius R — 2.0£ and k = 1.0 for 
— £/& = 0.0 (solid curves), 0.2 (dashed curves) and 0.4 
(dash-dotted curves). The thick curves indicate the supercon- 
ducting/normal transition and the thinner curves the ground 
state transitions between the giant vortex states with different 
vorticity L. 



FIG. 1. (a) The ground free energy and (b) the magneti- 
zation as a function of the applied magnetic field for a su- 
perconducting cylinder with radius R = 2.0£ for k = 0.28 
and — £/& = (solid curves), 0.2 (dashed curves) and 0.4 
(dash-dotted curves). 



FIG. 10. The (—£/&) — H phase diagram for a supercon- 
ducting cylinder with radius R = 4.0£ for ft = 0.28 (a) and 
k — 1.0 (b). The solid curves indicate the transitions between 
the different Instates and the thick curve gives the supercon- 
ducting/normal transition. 



8 



FIG. 11. (a) Free energy and (b) magnetization for vortex 
states in a superconducting cylinder with radius R = 4.0£ and 
k = 1.0 for the normal boundary condition — £/& = 0. Giant 
vortex states are given by solid curves, multivortex states by 
dashed curves and the open circles indicate the multivortex 
to giant vortex transition fields. The vertical lines in (b) give 
the ground state transitions. 

FIG. 12. (a) Free energy and (b) magnetization for vortex 
states in a superconducting cylinder with radius R — 4.0£ and 
k = 1.0 for the normal boundary condition — £/& = 0.2. Giant 
vortex states are given by solid curves, multivortex states by 
dashed curves and the open circles indicate the multivortex 
to giant vortex transition fields. The vertical lines in (b) give 
the ground state transitions. 

FIG. 13. Free energy for the L = 3 state in a superconduct- 
ing cylinder with radius R = 4.0£ and k = 1.0 for — £/& = 0.0, 
0.2, 0.4, 0.6, 0.8 and 1.0. Giant vortex states are given by solid 
curves, multivortex states by dashed curves and the open cir- 
cles indicate the multivortex to giant vortex transition fields. 

FIG. 14. The Cooper-pair density for the multivortex state 
with L = 3 in a cylinder with R/£ = 4.0 and n = 1.0 at 
H /H c2 = 0.445 for -£/& = 0.0 (a), 0.2 (b), 0.4 (c), 0.6 (d), 
and the logarithm of the Cooper-pair density for — £/& = 0.8 
(e) and 1.0 (f). 
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